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ABSTRACT 

Astrometric monitoring of the nearby early-L dwarf DE0823-49 has revealed a low-mass companion in a 248-day orbit that was 
announced in an earlier work. Here, we present new astrometric and spectroscopic observations that allow us to characterise the system 
in detail. The optical spectrum shows Li i-absorption indicative of a young age and/or substellar mass for the primary component. The 
near-infrared spectrum is best reproduced by a binary system of brown dwarfs with spectral types of LI.5 -t L5.5 and effective 
temperatures of 2150 ± 100 K and 1670 ± 140 K. To conform with the photocentric orbit size measured with astrometry and the 
current understanding of substellar evolution, the system must have an age in the 80-500 Myr range. Evolutionary models predict 
component masses in the ranges of Mi =: 0.028 - 0.063 Mq and M 2 - 0.018 - 0.045 Mq with a mass ratio of g 0.64 - 0.74. Multi¬ 
epoch radial velocity measurements unambiguously establish the three-dimensional orbit of the system and allow us to investigate its 
kinematic properties. DE0823-49 emerges as a rare example of a nearby brown dwarf binary with orbit, component properties, and 
age that are characterised well. It is a juvenile resident of the solar neighbourhood, but does not appear to belong to a known young 
association or moving group. 
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1. Introduction 

Binary stars are vital for advancing our understanding of stel¬ 
lar formation and evolution. They provide us with the oppor¬ 
tunity of direct mass measurements through orbital motion, and 
they are calibration systems for evolutionary models, which have 
to match the observations of a two-body system that presum¬ 
ably is co-eval and has near-identical composition. Very low- 
mass binary systems composed of brown dwarfs and/or ultracool 
dwarfs (spectral type M7 and later) are no exception and are in 
many ways more enlightening, given that brown dwarfs dim and 
cool substantially over time, yet they are less common than their 
higher-mass counterparts. The discovery and characterisation of 
ultracool dwarf binaries is thus a rare opportunity to yield obser¬ 
vational input for the refinement of models, eventually leading 
to a better understanding of ultracool dwarf physics. 

The L dwarf DENIS J082303.1-491201, hereafter 
DE0823-49, was targeted as part of an astrometric planet 
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090.C-0786, and 092.C-0202. 
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search (|Sahlmann et al.|2014|l and consequently identified 


as an 


astrometric binary (Sahlmann et al. 


2013|l. On the basis of fea¬ 


tures in its optical spectrum ( |Phan-Bao et al.||2008| l, |Sahlmann| 
|et al.| ( |2013| l suggested that DE0823-49 is younger than the 
average field population, which corresponds to a primary at 
the hydrogen-burning mass limit orbited by a low-mass brown 
dwarf secondary. Here, we present follow-up spectroscopic 
observations that allow us to better constrain the system’s 
properties, in particular the components’ effective temperatures, 
the system age, and individual masses. 


2. Observations 

2.1. VLT/FORS2 imaging 

We obtained two new astrometric epochs with the EORS2 instru¬ 
ment ( [Appenzeller et al!]|1998| l at ESO’s Very Large Telescope 
(VET), in addition to the data used for the orbit discovery. De¬ 
tails on the observation strategy can be found in |Sahlmann et al.| 
P013||20T4| |. The photocentre computations forDE0823-49 are 
complicated by a close background star, which resulted in no¬ 
ticeable systematic errors for this object compared to other tar¬ 
gets of the astrometric survey. Therefore a special model that 
extends the work of [Lazorenko et al.| ( |2014| l was devised for 
the reduction of DE0823-49 images. In this model, the seeing- 
dependent light contribution of the background star is taken into 
account. A second improvement concerns the enhanced suppres¬ 
sion of parameter correlations when reducing data of objects that 
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show significant orbital motion. All data were re-reduced with 
the improved methods and used for the orbit adjustment. 


2.2. VLT/UVES spectroscopy 

We observed DE0823-49 on 2013 October 7 (MJEp 
56572 .3332) with the red arm of UVES at the VET ( |DekkeretaL 


2000 1 using a 172 slit width, which provided a resolving power 
of R~33 000, and the Dichroic 2 standard setup centred at 760 
nm to cover the wavelength range of 565 - 931 nm. The expo¬ 
sure time was 2830 s and the observation took place with I'.'Ol 
optical seeing at an airmass of 1.52. The spectrum was recorded 
on two separate chips (REDE and REDU) and was reduced using 
the ESO pipeline in standard setup. 


2.3. IRTF/SpeX spectroscopy 

We observed DE0823-49 on 2013 November 24 (UT) with the 
SpeX spectrograph on the NASA Infrared Telescope Eacility. 
(IRTE; |Rayner et al.||2003) . Conditions were clear but windy 
with poor seeing (174 at //-band). We used the SpeX prism mode 
with the 075 slit aligned with the parallactic angle, yielding 0.8- 
2.45 pm spectra with an average resolution A/AA « 120. Eight 
exposures of 120 s each were obtained at an airmass of 2.81, fol¬ 
lowed by observations of the AOV star HD 83719 (V = 7.64) 
at an airmass of 2.79. HeNeAr arc lamps and quartz lamp expo¬ 
sures were also obtained for dispersion and pixel response cali¬ 
bration. Data were reduced using SpeXtool version 3.4 ( | Cushing] 
jet al.|2004{ jVacca et al.|2003| ) following standard procedures for 
point-source extraction. 


2.4. Keck/NIRSPEC spectroscopy 


We obtained two high-resolution infrared spectra of DE0823-49 
with the NIRSPEC echelle spectrograph on the Keck II tele¬ 
scope (McLean et al. 2000| l on 2014 April 14 and December 
8. Conditions on both nights were clear with 076 and 079 see¬ 
ing at X-band, respectively. We used the N7 order-sorting filter 
and 07432-wide slit to obtain 2.00-2.39 pm spectra over orders 
32-38 with A/AA = 20000 (Av = 15 km s“^) and dispersion of 
0.315 A pixel '. Two dithered exposures of 600 s each were ob¬ 
tained at an airmass of 2.78, followed with observations of the 
AO V star HD 87363 (V = 6.11). Elat held and dark frames were 
obtained at the start of the night with the same instrument setting. 


3. Analysis and interpretation 

Because of the new astrometric measurements and the improved 
reduction method, we re-analysed the astrometric data that com¬ 
prise 16 epochs over a timespan of 831 days. The analysis meth¬ 
ods are described in detail in jSahlmann et al.| ( |2013| l and consist 
of a genetic algorithm followed by a Markov-Chain Monte Carlo 
(MCMC) code. The updated ht parameters and their conhdence 
intervals are reported in Table where AaJ and A5q are relative 
offsets to the target’s position at the reference date TRef taken as 
the arithmetic mean of the observation dates, m is the relative 
parallax, pa* and ps are the proper motions, e is the eccentricity, 
to is the argument of periastron, P is the orbital period, dRef is 
the mean longitude at TRef, Q is the ascending node, i is the or¬ 
bit’s inclination, and a is the semi-major axis of the photocentric 
orbit. The parameters p and d model the differential chromatic 

* Modified Julian date (MJD) is barycentric Julian date - 2400000.5. 
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Table 1. Updated astrometric parameters of the DE0823-49 system. 


Atto 

(mas) 

-235.04+“|^ 

Ado 

(mas) 

-13.19!«4 

w 

(mas) 

48.27!°;j2 

Pa> 

(mas yr"^) 

— 1 Q9+0.06 

Pd 

(mas yr“') 

7 QQ+o.oe 
' '^^-0.06 

e 


^■■^^-0.04 

CO 

(deg) 

41 8+"^'* 

p 

(day) 

947 7S+0-64 
j_0.64 

dRef 

(deg) 

-365.8!20 

fi 

(deg) 

-13.8 + 2.0 

i 

(deg) 

52 2+'-"^ 

a 

(mas) 

4.62!°;[2 

p 

(mas) 

20.8![;^ 

d 

(mas) 

-26.2+ 1.5 

Sq, 

(mas) 

0 15+^-^^ 

^ '^-O.lO 

Sd 

(mas) 

^■^^-0.08 

Derived and additional parameters 

PRef 

(MJD) 

55926.823928 

Am 

(mas) 

-0.06 + 0.04 

^abs 

(mas) 

48.33 + 0.14 

Distance 

(pc) 

20.69 + 0.06 

Number of epochs / frames 

16/334 

^ 0-C,Epoch 

(mas) 

0.176 


Notes. Parameter values are the median of the marginal parameter dis¬ 
tributions and uncertainties represent Icr-equivalent ranges. The proper 
motions are not absolute and were measured relative to the local refer¬ 
ence frame. The inclination is measured relative to the sky plane. 


refraction and Sa and sg are nuisance parameters. The parallax 
correction Am was determined in jSahlmann et al.j ( |2013| ) and 
yields the absolute parallax lUabs- The time of periastron passage 
To can be retrieved via the mean anomaly 

M = /l-w = 27r(^-,^o) ^ To = TRef-P^, (1) 

where t is time and 0o = To/P is the phase at periastron. 

The most significant change is the reduction of the resid¬ 
ual scatter from 0.330 milli-arcseconds (mas) ( jSahlmann et al.j 
2013|l to CTo-c.Epoch - 0.176 mas, which reflects the improved 
data reduction. Consequently, the confidence intervals of most 
parameters are smaller as well. A notable change in the orbital 
parameters is a smaller inclination (56.6+2.0° before, 52.2+1.5° 
in this work). Other parameters, in particular the parallax, pe¬ 
riod, eccentricity, and photocentric semi-major axis a, remain 
unchanged when accounting for the uncertainties. 


3.1. Radial velocity measurements 


The Keck/NIRSPEC spectra were optimally extracted using a 
modified version of the REDSPEC package; the spectra on each 
night had a S/N~12. These data were then forward-modelled us¬ 
ing a custom MCMC implementation of the method described 
Blake et al.|( |2010| l. We used the Solar atlas of jLivingston &] 


Wallace (19911 to model telluric absorption and the BT-Settl at¬ 
mosphere models ( Allard et al.|2011) to model the spectrum of 
DE0823-49. A T^ff = 2000 K, logg = 5.0 (cgs) model provided 
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Offset in Right Ascension (mas) 

Fig. 1. Astrometric motion of DE0823-49 and its photocentric orbit, 
updated from |Sahlmann et al.|p013y Panel a shows proper and paral¬ 
lactic motion relative to the field of reference stars. Panel b is a close-up 
of the photocentric orbit caused by the gravitational pull of the orbiting 
brown dwarf. Observations with uncertainties and the best-fit model are 
shown as black circles and grey curve, respectively. 


the best fit. Figure]^ shows that the extracted spectrum is a good 
fit to the final model. The distribution of chain values yields the 
mean heliocentric radial velocities and projected rotational ve¬ 
locities listed in Tabled 


Wavelength (pm) 


2.29600 2.29957 2.30306 2.30646 2.30979 2.31304 2.31622 2.31934 



Pixels 

Fig. 2. High-resolution K-band spectrum of DE0823-49 obtained with 
NIRSPEC (black line) compared to a best-fit spectral (red line) and 
spectral plus telluric model (green line). Dilference between data and 
model is plotted in black at bottom; the uncertainty spectrum is indi¬ 
cated in grey. The values listed below the target name on the top right 
are Teff/log g/RV/v sin i. 

We also measured the radial velocity on the UVES spectrum 
using four strong atomic lines clearly visible in the data, namely 
Rbi at 780.0268 nm, Rbi at 794.7603 nm, Csi at 852.1132 nm, 
and Cs i at 894.3474 nm. The rest wavelengths for these lines 
were obtained from the NIST Atomic Spectra Database (version 


5.1; |Kramida et ar||2012[ ). The wavelengths of the centroid of 
the lines in the observed spectrum were derived using line pro¬ 
file fitting with the task splot in IRAF. Heliocentric correction 
was applied using the IRAF task rvcorrect and the information 
of Julian Date provided in the FITS header. This procedure was 
checked using archival UVES spectra for the FI dwarf 2MASS 
J1048428 1 -hOI 11580 from ESO program ID 078.C-0025(A) (PI 
A. Reiners). With our method we obtained a heliocentric radial 
velocity of 21.9 + 1.4 km/s, which is consistent with the pub¬ 
lished value of 24.0 +1.1 km/s ( jSeifahrt et al.|p0T0l l. No sig¬ 
nificant (i.e. larger than our uncertainties) radial velocity offset 
is expected to exist between the NIRSPEC and UVES measure¬ 
ments, as these two instruments have been shown to give compa¬ 
rable results for the brown dwarf LP944-20, although with larger 
scatter in the UVES based measurements ( |Martm et al.||200^ . 
For DE0823-49, we determined a heliocentric radial velocity of 
+6.5 ±1.5 km/s. 

A fourth radial velocity of +12.5 + 2.4 km/s was obtained 
from observations with the Magellan Echellette Spectrograph 
(MagE) in 2009. This measurement will be discussed in Bur- 
gasser et al. (in prep.). Using 32 sources that were both observed 


with MagE and with UVES (Seifahrt et al. 


2010 Reiners & 


|Basri|2009^ we determined an average velocity offset between 
both instruments and reduction methods of -0.29 + 0.59 km/s. 
It is small and not significant, especially when compared to the 
uncertainty of the MagE measurement of DE0823-49, and we 
therefore neglected it in our analysis. 


3.2. Spectroscopic age indicators 

Figure|^shows the reduced SpeX spectrum of DE0823-49 com¬ 
pared to four sources with reasonable age estimates and similar 
spectral types that were chosen to best match the 1.3 - l.Apm 
water band; TWA 29 (~8 Myr, i.e. the age of the TW Hy- 
drae Association), 2MASS J06085283-2753583 (2M0608-27, 
~20 Myr, i.e. the age of the Beta Pic Moving Group), EROS- 
MP J0032-4405 (EROS 0032-44, -100 Myr, i.e. the age of 
the Pleiades Moving Group or local association) and DENIS 
J170548.3-051645 (DE1705-05, -500 Myr, i.e. the age of the 
Ursa Majoris Moving Group). The comparison sources are from 
[Allers & Liu| ( |201^ and the references therein. 

DE0823-49 does not look like a source with the young age 
of TWA or Beta Pic, whose members lack FeH features around 
1 pm and have different H and K band shapes. The best match 
seems to be EROS 0032-44, i.e. an age of perhaps -125 Myr 
(the age of the Pleiades Moving Group is still debated, see 
Famaey et al. 2008|l, a source which exhibits low-gravity fea¬ 


tures and lithium absorption (|Martm et al.|1999| l. This is consis¬ 
tent with the [Allers & Liu| ( |2013| l gravity-sensitive indices that 
point towards an intermediate gravity (INT-G), indicating an age 
estimate of roughly 100-300 Myr. The MagE spectra reported 
in Burgasser et al. (in prep.) exhibits subtle indications of low 
surface gravity — slightly weaker 8183/8195 A Nat lines and a 
modestly stronger 8100 A VO band than equivalently-classified 
dwarfs — but these are not particularly strong features and thus 
consistent with a source in the few 100 Myr age range. 

The UVES spectrum of DE0823-49 shows a strong Li i ab¬ 
sorption feature. A Lit equivalent width of 2.4 + O.IA was 
measured by direct integration of the line profile shown in 
Fig . 1^ Such a strong Li i line has already been reported in high- 
resolution spectra of late-M and L dwarfs, and it indicates that 
lithium has not been depleted in the primary component of the 
binary ([Pavlenko et al.|[2007|). The spectral resolution of our 
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Wavelength (|j.m) 


Fig. 3. SpeX spectrum of DE0823-49 (black) compared to four sources 
(red curves) with similar spectral types and well-constrained ages. 
Source identifier, optical spectral type, and the approximate age are 
indicated next to every spectrum. Cluster memberships (UMa = Ursa 
Majoris, PMG = Pleiades Moving Group, BPMG = Beta Pic Moving 
Group, TWA = TW Hydrae Association) are from |Allers & Liu|p013| l. 


UVES spectrum is much higher than the resolution that has been 
used to classify L dwarfs using gravity sensitive features ( |Cruz| 
et al.|2009l, and we refrain from putting our data in a classifica¬ 


tion system based on low-resolution spectra. In future work we 
plan to compare our UVES spectra of the |Sahlmann et al^p014) l 
targets with high-resolution synthetic spectra to obtain quantita¬ 
tive estimates of their surface gravities (Martin et al. in prep.). 

As discussed in Magazzu et al. ( 1993) 1, the presence of 
a strong Lii line coupled with the effective temperature of 
DE0823-49 already indicates that it must have a substellar mass 
lower than about 0.065 solar masses and a relatively young age 
between 100 and 1000 Myr. A more precise estimate on the age 
of DE0823-49 using lithium and other age indicators is pre¬ 
sented in Sect. 13.41 


3.3. Spectral binarity 


Unresolved ultracool binaries with late-M/L dwarf and T dwarf 
components can be discovered and characterised by disentan¬ 
gling the spectral features of individual components in their 
combined-light spectra (e.g. |Burgasser|2007]l. For DE08 23-49, 
the estimated component masses from Sahlmann et al. ( 2013| l 
for an age of ~1 Gyr suggested that the secondary of this system 
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Wavelength (nm) 


Fig. 4. The Lii line in the UVES spectrum of DE0823-49. The dotted 
line marks the integration limits and continuum level chosen to measure 
the equivalent width. 


could be a T dwarf. However, we do not see the //-band ‘dip’ 
feature commonly present in combined-light L/T binary spectra, 
only a slight excess in the 2.05 and 2.2 pm region. 

We nevertheless applied the methods described in |Burgas^er 
|et al.| ( |2OT0| ) and compared the spectrum of DE0823-49 to bi¬ 
nary templates constructed from L and T d warf spectra in the 
SpeX Prism Libraiyp](SPL, [Burgasser 2014 1 . This procedure al¬ 
lows us to estimate the spectral types and effective temperatures 
of the individual components and their /-band magnitude differ¬ 
ence. 

Because of the Li i absorption and its near-infrared spectrum, 
we know that DE0823-49 is relatively young. To reflect this 
in spectral binary fitting, we performed the analysis with three 
sets of templates: those with only ‘young’ sources, only ‘not 
young’ sources, and with ‘all’ sources. In each case, this anal¬ 
ysis uses binary templates constructed from M9-L5 primaries 
and L4-T6 secondaries in the SPL, where individual component 
spectra were scaled to the absolute T-band magnitude - spectral 
type relation of |Dupuy & Liu| ( f2012 i. The F test was applied to 
compare the quality of different fits, which relies on the assump¬ 
tion of Gaussian uncertainties The young sources were selected 
as having previously identified signatures of low surface gravity 
(e.g. Cruz et al.|2009| Allers & Liu|2013 l, however this sample 
is dominated by sources older than ~100 Myr. The ‘not young’ 
sample excludes these young sources. 

For both the ‘young’ and ‘all’ sets of templates, the best sin¬ 
gle match to the spectrum of DE0823-49 is that of the LI.5 
2MASS J20575409-0252302 (hereafter 2M2057-02;|Cruz etal. 


2003[ Burgasser et al.|2004] Allers &Liu|2013| l, which is shown 
in Fig. 1^ The best-fit binary template for the unrestricted tem¬ 
plate set (‘all’) is shown in the top panel of Fig. Is i 


IS a com- 


^ http: //WWW. browndwarfs. org/spexprism __ 

^ Note that the used^“ prescription (Eq. (1) in Burgasser et al. |2010f 


does not account for the noise spectrum of the template, because this 
would bias the result towards the most noisy templates. 
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Fig. 5. Result of spectrum fitting with single templates. The single tem¬ 
plate 2M2057-02 (blue line) is the best fit to the SpeX spectrum of 
DE0823-49 (hlack line) both for the ‘young’ and ‘all’ set of templates. 
At the top of the panel, the spectral ranges considered for the fit are 
indicated by grey horizontal bars, where telluric water absorption dom¬ 
inates the gaps. The inset shows a close-up of the //-band region. The 
residual spectrum is shown in grey. 


bination of the L1.5 2MASSW J1411175-H393636 (Kirkpatrick 
et al.|2000|l and the L2.5 pec 2MASS J00145575-4844171 ( |Kirk ' 


Patrick et al.|2008 i, where the latter’s near-infrared appearance 


classifies it closer to a spectral type of L4. This provides a sig 
nificantly better match than the spectrum of 2M2057-02, par¬ 
ticularly in the shape of the 2.1 jjm K-band peak (likely due to 
stronger H 2 absorption in the secondary), so that the single tem¬ 
plate model is rejected with a confidence of 99.9 % on the basis 
of the F test. The influence of the companion on the blended light 
spectrum is thus clearly detected. 

For the ‘young’ template set, the bottom panel of Fig. 
shows the best-fit young binary template, a combination of 
DENIS-P J170548.4-051645.7 (DE1705-05, |Kendall et al. 
and 2MASS J20025073-0521524 (2M2002-05, |Cruz 


et al. 


|2007|l. The addition of the L6 fills in the excess flux at 


/T-band, and this binary template provides a significantly better 
match than the spectrum of 2M2057-02. The single template 
model is rejected with a confidence of 75 % on the basis of the 
E test. 

To obtain the parameters of DE0823-49 and their confidence 
intervals, we combine all template fits and assign them relative 
weights according to their values to compute weighted aver¬ 
age values with uncertainties like in |Burgasser et al. ( 2010| l. The 
results are summarised in Tablej^ which also lists the number of 
available binary templates and the E-test confidence. Effective 


temperatures were derived using the relations of Looper et al. 
( 2008| l and include an additional 0.5 spectral subtype uncertainty 
for optically classified sources. We measured the 2MASS J mag¬ 
nitudes for both components of the binary templates and used the 
[Schmidt et al.| ( 20T0] l relation between spectral type and SDSS i 
- 2MASS J colour to obtain estimates for the SDSS i magnitude 
difference between primary and secondary. We found that this 
is indistinguishable from integrating the Bessel I or the SDSS i 
bandpasses from the template spectra directly. The uncertainties 
are dominated by scatter in the viable sources, thus A/sdss - A/. 
Eor every individual spectral fit, we also computed bolometric 



Wavelength (pm) 


Fig. 6. Results of spectrum fitting with binary templates (cf. Fig. [^. 
Top'. When using ‘all’ templates, the SpeX spectrum of DE0823-49 
(black line) is best fit by the binary template (brown line) that is a 
combination of 2MASSW J1411175-t393636 (red line) and 2MASS 
J00145575-4844171 (blue line). Bottom'. The same when using only 
‘young’ templates. The best-fit binary template (brown) is a combina¬ 
tion of DE1705-05 (red) and 2M2002-05 (blue). 


corrections (Liu et al. 2010| l for the components of the spec¬ 
tral templates and estimated the /T-band magnitude difference 
AKmko- We used these to compute the component bolometric 
luminosities Lboi based on the system’s combined light absolute 
MKO magnitude, obtained from absolute 2MASS Kj and the 
MKO-2MASS correction from the spectrum. The uncertainty 
in Lboi includes 0.08 mag uncertainty in the /T-band bolometric 
correction. 


In all cases, the spectral types of the primary and secondary 
component are consistent with LI.5 and L5, respectively. All 
other properties are also derived consistently within their un¬ 
certainties for different template sets. Eor the ‘not young’ and 
‘air sets, the single template model is rejected with a confidence 
of 99.9 %, which we interpret as a significant detection of the 
companion in the blended light spectrum. Eor the ‘young’ tem¬ 
plates, the rejection confidence is 75 %, which would be consid¬ 
ered insufficient evidence of (young) multiplicity based on spec- 
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Table 2. Results of spectral binary fitting with different sets of templates 



‘young’ 

‘not young’ 

‘all’ 

# Primaries 

59 

381 

440 

# Secondaries 

28 

341 

369 

# Binaries 

1634 

128 868 

161121 

Primary SpT 

Ll.5 + 0.6 

Ll.l ±0.6 

Ll.l ±0.6 

Secondary SpT 

L5.5 ± 1.1 

L4.2+ 1.4 

L5.3 ± 2.3 

Primary Teff (K) 

2150+100 

2190 + 90 

2200 ± 90 

Secondary Tcir (K) 

1670+ 140 

1810+180 

1720 ± 230 

Primary logLboi (Eg) 

-3.74 ± 0.04 

-3.77 ± 0.07 

-3.73 + 0.08 

Second. logLboi (Lg) 

-4.20 ± 0.08 

-4.15 + 0.14 

-4.25 ± 0.20 

A/ (mag) 

2.4 ±0.6 

1.7 ±0.9 

2.3 ± 1.2 

F-test conf. (%) 

75 

99.9 

99.9 



troscopy alone. However, this statistic is skewed by the small 
number of young templates available in the SPL, and spectral 
binary ‘significance’ is less relevant since we already know that 
this system is a binary. We can still use the F-statistic to weight 
the spectral types from this analysis, which are in agreement with 
the ‘not young’ and ‘all’ templates samples. The high signifi¬ 
cance of the F test with ‘all’ templates may suggest that spectral 
binary fitting can be a viable procedure for detecting L-hL spec¬ 
tral type binary systems. We caution, however, that the binarity 
of DE0823-49 was established with astrometry, which simpli¬ 
fies the interpretation of the spectral fitting results. 

Because we know that DE0823-49 is relatively young, we 
adopt the parameters determined with the young templates. This 
yields component spectral types of Ll.5+0.6 and L5.5+1.1 and 
effective temperatures of 2150 + 100 K and 1670 + 140 K for 
the primary and secondary components of DE0823-49, respec¬ 
tively. The corresponding magnitude difference of AI = 2.4+ 0.6 
is valid for the passband of our EORS2 observations. 

3.4. Constraints on age and physical parameters 

The determination of individual effective temperatures in the 
previous section allows us to employ models of (sub-)stellar evo¬ 
lution to set constraints on the system’s age and its properties. 
The first comes from the detection of Li i absorption in the spec¬ 
trum, which implies a mass of < 0.065 Mq for the primary. Eig- 
ure shows how this translates into an upper age limit of ~0.5 
Gyr when coupled to the primary’s effective temperature and the 
DUSTY ( [Chabrier et al.|2000| evolutionary models. 

Additional constraints can be derived from the photocen¬ 
tric orbit, which was measured in the /-band (I_Bessel filter of 
EORS2, centred at 760 nm). Because of the moderate magni¬ 
tude difference between the two components of DE0823-49, we 
have to take into account the light contribution of the companion 
if we want to determine the barycentric orbit size that relates to 
the system’s mass function. The fractional mass 

/ - M2/(Mi + M 2 ) (2) 

and the fractional luminosity in the observation passband 

// = L2/(Ti+T2) = (1 + 10‘’'4^'")-', (3) 

where Am is the magnitude difference, define the relationship 
between the semimajor axis a of the photocentre orbit and the 
semimajor axis of the relative orbit, both measured in mas: 

a ^ Qrelif - P). (4) 

An independent constraint on the relative semimajor axis is 
given by Kepler’s law: 

G(Mi+M2)-4;r2^, (5) 
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Fig. 7. Upper age limit from the presence of Li i absorption in the pri¬ 
mary’s spectrum. The curve shows the effective temperature of objects 
that have depleted 90 % of their lithium (i.e. lithium abundance = 0.1) 
as a function of age after interpolation of the DUSTY models. Hori¬ 
zontal dashed lines indicate the primary’s effective temperature range, 
which translates into an upper age limit of ~0.5 Gyr. 



Fig. 8. Age constraints derived on the basis of two sets of evolutionary 
models, derived using effective temperature - mass relationships. The 
curves shows the photocentric orbit size ac for various ages and fixed 
effective temperatures, where the solid and dashed curve corresponds 
to the DUSTY and Saumon & Marley|p008| l models, respectively. The 
shaded regions delineate the uncertairities (dark grey for DUSTY, light 
grey for |Saumon & MarleyH2008) ). The horizontal lines show the mea¬ 
sured value of a and its uncertainty. Vertical dotted lines indicate age 
limits derived from the photocentric orbit. 


were G is the gravitational constant, is measured in metres 
and F is in seconds. The relation between and arei is given by 
the parallax. Any combination of the unknown parameters Mi, 
M 2 , and Am/ has to satisfy Eqs.|^and|^ which are constrained 
by the measured quantities F, a, and the parallax m. We used the 
DUSTY and |Saumon & Marley| ( [2008] ) models to impose theo¬ 
retical relationships between effective temperature or bolometric 
luminosity and mass as a function of age, which is the parameter 
we want to constrain. 

For every accepted binary template in the previous section, 
we derived theoretical component masses on the basis of the ef¬ 
fective temperatures and combined them with the estimated /- 
band magnitude difference and Eqs.[2j|^to compute the expected 
photocentric orbit sizes for a range of ages. For every age, we de¬ 
termined the effective orbit size ac and its uncertainty from the 
weighted average and the weighted standard deviation, respec¬ 
tively, where again we used the weights corresponding to the 
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Table 3. Derived age, masses, and mass ratio of DE0823-49 as a function of model and the relationship used to obtain theoretical masses, indicated 
in the column R. There are two rows for every configuration, which give the parameter ranges as delineated by vertical lines in Eigs.[^and[^and 
the best-fit values with uncertainties that reflect these ranges. 


Model 

R 

Age 

Ml 

M 2 

q — M 2 /M 1 



(Myr) 

(Mo) 

(Mo) 


DUSTY 

^bol 

[91,473] 

[0.028,0.060] 

[0.018,0.044] 

[0.65,0.74] 

DUSTY 

^bol 

216-^^^ 

^^"-124 

u.UHi_ooi3 


^-0.06 

SMOS" 

^bol 

[110,568] 

[0.031,0.063] 

[0.020,0.045] 

[0.65,0.71] 

SM08 

^bol 

041 +326 

-0.013 

u.uj)u_ooio 

0.68/°;“ 

SM08 

Tboi* 


u.UHZ,_ooi3 

^■^^^_0.009 

0.69/°;“ 

DUSTY 

Left 

[83,342] 

[0.034,0.054] 

[0.022,0.039] 

[0.64,0.73] 

DUSTY 

Left 

1 00+203 

u.u^i _0 006 

^■^^^_0.007 

0 70+003 
^-0.06 

SM08 

Left 

[86,428] 

[0.035,0.060] 

[0.022,0.043] 

[0.64,0.71] 

SM08 

Teff 

178/250 

0 045 
u.uhJq.oio 

U.UJ)U 0008 

0-681o:S 


Notes. 

a brown dwarf of a given mass 0.3 dex brighter in terms of logLtoi. 


Saumon & Marley 1 2008 1 . These results were obtained when considering a systematic offset in the model that corresponds to making 



Age (Myr) 

Fig. 9. Same as Eig. but derived using the bolometric luminosity - 
mass relationships. 


values of the spectral binary fits. When needed, the models were 
interpolated linearly in effective temperature, mass, or age. 

Figurej^shows the estimated orbit size ac a function of age. 
Using the DUSTY models results in slightly larger values com¬ 
pared to the |Saumon & Marley ( 2008| l models. The measured 
orbit size a yields an additional constraint on the system age. To 
address the effects of potential systematic differences between 
spectral type - effective temperate relations for young and nor¬ 
mal sources, which could be relevant in the case of DE0823-49, 
we repeated this analysis using the bolometric luminosities to 
obtain theoretical masses, leading to similar results, see Figure 

|9l 

Table lists the numerical results for all configurations, 
which give an overall consistent picture of DE0823-49 with an 
age of ~ 100-500 Myr, primary mass of Mi = 0.028 - 0.063 Mq, 
companion mass of M 2 = 0.018 - 0.045 Mq, and a mass ratio 
q — M 2 /Ml = 0.66 - 0.74. On average, the bolometric luminos¬ 
ity method yields a slightly older age range, but all results agree 
well with the upper age limit of ~500 Gyr set by the Li i detec¬ 
tion. A typical example of best-fit values is an age of 240 /j3q 
M yr, masses of Mi = 0.044/9 g[? and M 2 = 0.030/°°}^, and 


a mass ratio of q 


= 0 . 68 -°;“ 


(obtained using 


Saumon & Mar¬ 


ley ( 2008|l and bolometric luminosities.) These age estimates are 


consistent with the INT-G surface gravity classification and the 


quality of gravity-sensitive features present in the near-infrared 
and optical spectra. 

Einally, the age and mass parameters that we derived for 
DE0823-49 fall into a range where the results of |Dupuy et al.| 
( 2014) 1 suggest that theoretical models underpredict the lumi¬ 
nosities of brown dwarfs. To quantify the effect that such sys¬ 
tematic errors would have on our results, we introduced an 0.3 


dex offset in the Saumon & Marley (2008 1 bolometric luminosi¬ 
ties that would make a BD with given mass and age ss2 times 
brighter in LboU an amplitude that corresponds to observations 
( |Dupuy et al.|2009} . When including the offset, the resulting ac¬ 
ceptable age range of DE0823-49 is shifted to 150-860 Myr, but 
the component masses remain the same, see Table Because an 
upper age limit is set by the Li i detection, we conclude that such 
a systematic model error has negligible effect on the age and 
mass properties we derive for DE0823-49. 


3.5. Rotation 

Using the NIRSPEC infrared spectra, we measured a mean pro¬ 
jected rotational velocity of vsinf = 32 + 4 kms^', which is 
comparable to that found for other field LI.5 dwarfs ( {Reiners | 
& Basri||2008[). If the binary orbit orientation (sin/= 0.79) co¬ 


incides with the spin axis of the primary, the actual rotational 
velocity is ~40 km s *. We performed a periodogram analysis of 
the optical variability of DE0823-49 shown in Eig. 10 but found 
no evidence for periodic variations at correspondingly short ro¬ 
tation periods (^3 h for a Jupiter-sized body). We only noticed a 
slight variation of ~10 mmag over 800 days. 


3.6. Radial velocity orbit 

Eigure [TT| shows the radial velocity measurements (see Table 
and the expected radial velocity curve of DE0823-49A shifted 
by the systemic velocity y = 2.98 + 0.75 kms^' that we mod¬ 
elled as a constant offset to the measurements. We assumed 
masses of Mi = 0.044 Mq and M 2 = 0.030 Mq, but the curve 
is essentially the same for a younger or older configuration, e.g. 
Ml = 0.031 Mq and M 2 = 0.020 Mq. The available measure¬ 
ments are compatible with the shown curve, which allows us to 
conclude that we resolved the ambiguity of 180° in inclination 
and longitude of ascending node, which means that the values 
reported in Table are accurate. More radial velocity measure- 
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Fig. 10. Differential magnitude variation of DE0823-49 in /-band as a 
function of time. Grey symbols correspond to measurements in individ¬ 
ual FORS2 frames, whereas black circles show the epoch average. 


ments would be needed to independently constrain the system 
parameters. 



Year 



Orbital phase 


Fig. 11. Estimated radial velocity curve of DE0823-49A (Mi = 
0.044 Mq, M 2 = 0.030 Mq) as a function of time (top panel) and orbital 
phase (bottom panel, phase 0 = periastron). The measurements with 
UVES (solid circle), NIRSPEC (triangles), and MagE (open circle) are 
shown and the systemic velocity is indicated by the dashed line. 


4. Discussion 

Young sources in the solar neighbourhood are often found in 
associations that have similar kinematic properties. Since all 
relevant quantities, i.e. coordinates, parallax, proper motions, 
and radial velocity, are now known for DE0823-49, we exam¬ 
ined its potential membership in a young association using the 
BANYAN II web tool ( |Gagne et al.|2014| l. According to BANYAN 


Table 4. Radial velocity measurements of DE0823-49A 


Instr. 

Date 

Epoch 

(MJD) 

RV 

(km/s) 

V sin i 
(km/s) 

MagE 

2008-01-08 

54838.21756 

12.5 + 2.4 


UVES 

2013-10-07 

56572.33321 

6.5 + 1.5 

26 + 5 

NIRSPEC 

2014-04-13 

56760.22524 

9.7+ 1.0 

34 + 5 

NIRSPEC 

2014-12-08 

56999.55533 

8.4+ 0.8 

29 + 6 


II, DE0823-49 has a 100 % probability of belonging to the field 
population. The probability of belonging to one of the 7 young 
associations probed is reported as 0.00 %. 

DE0823-49 thus appears isolated in terms of membership, 
although there may exist numerous very-low-mass few-member 
young associations that are not recognised yet in the solar neigh¬ 
bourhood. Searches for young very-low-mass stars outside star¬ 
forming regions indicate that they are typically near larger star¬ 
forming regions but not inside them ( [Martin & Kun|199^ jVal-j 
jdivielso et al.|2009] l. 

Another consequence of the new age estimate is related to 
the mass ratio (q - M 2 /M 1 ) of the DE0823-49 binary that had 
a reported value of 0.36 derived under the assumption of a sys¬ 
tem age of ~1 Gyr ( jSahlmann et al.|2013| l. With the character¬ 
isation presented here, the primary mass is lower but the mass 
ratio is higher in the range of ~0.66-0.74. This system appears 
now more similar to the bulk of known ultracool dwarf bina¬ 
ries, whose mass ratio distribution is dominated by systems with 
q > 0.7 ( Burgasser et al.||2007| l. However, its mass ratio is still 
uncommon and lies in the tail of the ^'-distribution of field brown 
dwarf binaries (e.g. jLiu et al.||2010 l. The estimated companion 
mass in this binary system remains in the range of ~30 Jupiter 
masses, which corresponds to the most massive planets found 
around Sun-like stars (jSahlmann et al.|201 Ijl. 


5. Conclusions 

We have presented the detailed characterisation of the tight 
brown dwarf binary system DE0823-49. After the orbit discov¬ 
ery with optical ground-based astrometry, we obtained follow-up 
optical and infrared spectroscopy. The combination of spectral 
modelling with age indicators and evolutionary models leads us 
to the following conclusions: 

- DE0823-49 is a system composed of two brown dwarfs with 
spectral types of El.5+0.6 and L5.5+1.1 that orbit each other 
in ~ 248 days. 

- The spectrum of this system is signihcantly better matched to 
a binary template than a single spectral template, providing 
an indirect detection of the companion in the near-infrared. 
Based on the classification of the templates, we infer the 
components’ effective temperatures to be 2150 ± 100 K and 
1670 ± 140 K. 

- We estimate the age of the system to be in the range of 80- 
500 Myr, which is supported by the detection of Li i absorp¬ 
tion in the optical spectrum and by the properties of the near- 
infrared spectrum. Rehned upper and lower age bounds were 
derived from the measured photocentric orbit size and theo¬ 
retical relationships between mass and effective temperature 
or bolometric luminosity. 

- DE0823-49 does not appear to belong to any of the known 
nearby young associations. 

- Evolutionary models predict component masses in the ranges 
of Ml 0.028 - 0.063 Mq and M 2 - 0.018 - 0.045 Mq with 
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a mass ratio of ^ = 0.64 - 0.74. Both components are thus 
substellar objects below the lithium-burning mass limit and 
may have masses comparable to the most massive planets 
found around Sun-like stars. 

At a distance of 20.7 pc, DE0823-49 is a rare example of a 
nearby brown dwarf binary with well-characterised orbit, com¬ 
ponent properties, and age. 
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